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Ground-State  Properties  of  a  He  Impurity  in  Liquid  He  Monolayers 


Chung- In  Um  and  Sahng-Kyoon  Yoo 
Department  of  Physics 
College  of  Science 
Korea  University 
Seoul  136-701,  Korea 

Thomas  F.  George 

Departments  of  Chemistry  and  Physics 
Washington  State  University 
Pullman,  Washington  99164-1046 


We  present  the  results  of  ground-state  energies,  radial  distribution  functions, 
liquid  structure  functions  and  effective  interactions  for  a  ^He  impurity  in  a  ^He 
background  in  two  dimensions.  The  hypernetted- chain  scheme  for  the  system  described 
by  a  Jas trow- type  wavefunction  is  used,  taking  into  account  the  triplet  correlations 
'nd  elementary  diagrams  up  to  fifth  order.  Solving  the  Euler- Lagrange  equations  for 
the  two-body  distribution  functions,  which  contain  triplet  correlation  and  elementary 

3 

diagrams,  improves  the  results  considerably.  Furthermore,  as  a  He  impurity  is 

4 

inserted  into  the  He  background,  the  ground-state  energy  increases,  but  the 

equilibrium  density  decreases  from  0.0350  A  to  0.0336  A  .  The  radial  distribution 

function  is  broadened,  while  its  maximum  is  lowered  and  shifted  to  the  right  (the 

direction  of  increasing  radial  distance)  due  to  its  larger  zero -point  energy,  with 

3 

therefore  less  localization  of  the  He  particle.  The  results  are  compared  with  Monte 


Carlo  results  and  other  studies. 


PACS  No.:  67.60. -g,  67.60.Fp,  67.40.Db 


2 


I.  Introduction 

In  the  past  three  decades  there  has  been  considerable  interest  in  single- 

1  2 

component  Bose-fluid  and  binary  boson  mixtures  at  very  low  temperatures,  such  as 

3 

stable  bulk  spin-polarized  hydrogen  atom  and  its  isotope.  Among  the  Bose  systems, 

4  3  4 

liquid  He  and  He-  He  mixtures  have  been  analyzed  successfully  by  introducing  a 
Jastrow-type  wavefunct ion .  In  the  investigation  of  ground  states,  the  Green's 

function  Monte  Carlo  (GFMC)  technique,^  "*  paired-phor.cn  analysis  ^  average 

8  .9-10 

correlation  approximation  (ACA)  and  hypernetted-chain  (HNC)  approximation  scheme 

provide  powerful  tools  in  both  two  and  three  dimensions.  First,  three  calculations 

give  the  exact  ground-  and  excited-state  energies,  but  less  information  about  the 

optimized  correlation  function  and  other  quantities  which  depend  on  the  long-range 

behavior,  although  much  experimental  data  about  them  have  been  reported.  However,  in 

spite  of  its  weakness  in  the  calculation  of  the  ground-state  energy,  the  HNC 

approximation  can  be  improved  by  including  the  triplet  correlation  function  and  some 

11  12-17 

elementary  diagrams  in  the  energy  and  Euler- Lagrange  equations.  Woo  and 

12-13 

Coldwell  have  considered  the  wave  function  consisting  of  the  Jastrow  function 

multiplied  by  a  three-particle  function  and  obtained  the  improved  energy  and  liquid 

structure  function.  Recently,  Fabrocini  and  Polls^  have  extended  this  idea  to  the 

ideal  boson  He-  He  mixture  in  the  zero  He  concentration  limit. 

18 

In  the  two-dimensional  case,  Miller  and  Woo  have  studied  the  ground  state  of  a 
^He  monolayer  assuming  no  significant  modifications  by  the  existence  of  a  substrate 
in  the  mobile  limit.  Chang^  has  evaluated  the  ground-state  energy  and  liquid 
structure  function  using  the  optimal  Jastrow  function  determined  from  a  self- 
consistent  PPA,  and  Hatzikonstantinou^  has  also  investigated  ground-state  properties 
by  using  a  Jastrow-type  wavefunction  and  solving  a  Euler -Lagrange  equation  within  the 


HNC  approximation  without  a  three-particle  factor. 


3 


The  main  purpose  of  the  present  paper  is  to  evaluate  the  ground- state  properties 
434  3 

of  two-dimensional  He  and  He-  He  mixtures  in  the  zero  He  concentration  limit  by 
optimizing  a  Jastrow-type  wavefunction  and  solving  the  Euler-Lagrange  equations 
within  the  HNC  approximation,  including  a  triplet  correlation  function  and  elementary 
diagrams  up  to  fifth  order. 

In  Sec.  II,  we  summarize  the  construction  of  the  Euler-Lagrange  equations  for 
3 

zero  He  concentration  obtained  from  the  energy  minimization  condition.  The  Jastrow 

wavefunction  consisting  of  the  product  of  the  two  and  three-body  correlation  factors 

and  the  Jackson  and  Feenberg  energy  form  are  adopted.  We  apply  our  results  to  the 
4  3  4 

He  monolayer  and  the  He-  He  mixture  in  Sec.  Ill  to  obtain  the  ground-state 
energies,  radial  distribution  functions,  liquid  structure  functions  and  the  effective 
potential,  assuming  a  Lennard-Jones  6-12  potential.  The  improvements  in  the  ground- 
state  energies  and  other  quantities  are  also  discussed.  In  Sec.  IV,  we  make 
conclusions . 


II.  Euler-Lagrange  Equation  In  the  Zero-Concentration  Limit 

In  this  section  we  consider  a  homogeneous  binary  boson  mixture,  which  can  be 
described  the  the  Hamiltonian 
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Here,  the  first  two  terns  are  kinetic  energies  of  components  a  and  /? ,  respectively, 
and  V^(r)  is  the  pairwise  interaction  potential  for  the  three  different  pairs  of 
particles.  We  assume  that  the  total  number  of  particles  in  an  area  A  is  N  -  Nq  +  N^. 
We  choose  the  Jastrow-type  wavefunction  in  terms  of  two-  and  three  -  particle 
correlation  functions  of  the  form 
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where  u_^(r)  is  the  two-body  correlation  function  which  satisfies  the  boundary 
conditions 


lira  u  (r)  -*■  «> 
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lim  UQfl(r)  “  0  ' 

r*+<o  ” 


(2.3) 


and  ua^(r)  is  c^e  three-body  correlation  factor.  The  two-body  distribution  function 
can  be  defined  as 
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where  dr(i  ,  i„)  denotes  dr,  .  ..dr„  with  the  exclusion  of  dr.  and  dr.  .  and  p  -  N  /A 
'  a  B  IN  l  l  _  ct  a 
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(2) 


is  the  partial  density  of  component  a.  The  radial  distribution  function  g  a  (r)  must 

Qp 

also  satisfy  the  conditions 


lim  [g  (r)  -  1]  -  0 
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P  I  d  r  Iga/J(r)  -  1]  +  Saf) 


(2.5) 


Defining  x  —  p  /p  as  the  concentration  of  component  a,  where  p  —  (N  +  N_)/A  is  the 
a  a  a  p 

total  density,  we  can  express  the  energy  per  particle  of  the  system  as 
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(3)  32  23 
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The  two-  and  three -body  energies  are 
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with  the  effective  interaction 
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where  g^^(r  ,rrt,r  )  is  the  three-body  distribution  function.  This  is  written  as 
apy  a  p  7  y 
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where  A(r  ,  r-,r  )  is  the  sum  of  the  contributions  of  all  the  Abe  terms. 

P  7 

(2)  (2) 

In  order  to  establish  a  relation  between  g^  (r)  and  f^  (r),  we  introduce  the 
HNC  approximation,  which  consists  of  the  equation^- 


-  [f^><r>]2exp[C,(S<r)  +  T^(r)  +  E^(D]  . 


(2.14) 


where 
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T  (r)  is  the  integral  of  the  triplet  function 
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(2.16) 


and  E  (r'l  i<?  the  sum  of  all  the  elementary  diagrams.  Rewriting  Eq.  (2.9)  and  taking 
q/3  ' 

C  a(r)  in  momentum  space,  we  get 
QP 


1  f  &  (2)/rW  ,r*  _B?A  f  d2r  B(2)f  n-2,  ,(2), 

2  P  J  d  r  &af}  (r)%(r)  '  8m  j  U  r  go/3  v~'*  ""ba/3 


M2p  So3 _ 

8ra  a  .2 

q/3  (2rr)  PqP^ 


J  o2*  <Vk)[Wk) 


S_,<k)l'  -  K »(Wi  K2 


1  -  8^  I  d2?  8")(r)V2|T^(r)  ♦  E  CD]  .  (2.17) 
a/3 


where  S  (k)  is  the  liquid  structure  function  of  the  form 
op 
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Vk)  "  ^ a/3  +  2np  J  dr  r  ^>(r)  -  L)  J0(kr)  •  (2‘18) 

2 

D(k)  —  ^h^22  *  ^12  and  Q  **  7-  Here  J^(kr)  is  the  zeroth-order  Bessel  function.  To 

(2) 

minimize  the  energy,  we  must  calculate  the  variational  equations  6E/6fv  (r)  -  0  and 
(3) 

SE/Sf  (r)  -  0  and  solve  the  resultant  Euler- Lagrange  equations.  In  the  presence  of 
triplet  correlations  this  is  very  complicated.  To  solve  this  problem  we  take  the 
limit  of  zero  concentration  of  component  a  and  assume  that  f'  ;(r)  is  fixed  for  all 
the  triplets.  Minimizing  the  energy,  we  finally  get  the  Euler- Lagrange  equations 
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(2.19) 


0  T  E 

where  a  -  1,2.  Here  W  „(r),  W  _(r)  and  W  „(r)  are  the  induced  potentials  due  to  two- 

aZ  aZ  a Z 

body  correlations,  triplet  correlations  and  elementary  diagrams,  respectively. 
( r )  is  expressed  in  momentum  space  as 
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Here  we  ass-aue  E^2(r)  "  E2^(r).  We  may  rewrite  Eq .  (2.6)  for  slight  changes  in  as 
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Assuming  that  there  is  an  impurity  in  the  mixture,  then  x^  would  be  1/N,  and 

(2) 

approximately  to  order  1/N  we  can  minimize  E  through  minimization  of  E  2(r)  ,  that 

(2) 

is  g  (r)  can  be  estimated  consecutively. 

From  the  Euler-Lagrange  equations  we  can  construct  the  n-th  step  iteration 

(2)  (2)  (2)  (2) 

scheme  and  find  g  .  ,,(r)  -  g  0  (r)  +  6  g  ^  (r)  for  a  given  g  (r).  Then  by 

QZ,n+i  az,n  a2,n  al  ,n 

means  of  the  convolution  theorem,  the  Euler-Lagrange  equations  can  be  converted  to 
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where  W*  is  given  by 
etc 


W  - 

12 


k2[2m1'2m2S22  +  (1 


2m12m2)S22  -  1] 


22 


W  - 

22 


2  3 

k  ( s 
K  -  22 


-  1) 


2S 


22 


(2.25) 


10 


This  variational  calculation  consists  of  two  parts:  (a)  we  consider  a  system  of  a 

background  component  and  a  trial  wavefunction  of  the  form  of  Eq .  (2.2)  and  determine 

(2) 

^22  variationally ,  and  (b)  we  replace  the  background  system  with  an  impurity  and 

introduce  the  trial  wavefunction 
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(2) 

Through  step  (b) ,  g^  (r)  can  be  determined  variationally . 

Equation  (2.24)  represents  a  set  of  integro-dif ferential  equations.  Under  the 

(2) 

assumption  that  Sg^  (O  vanishes  outside  some  radius  R,  Eq .  (2.24)  can  be  replaced 

by  a  finite  sum  over  a  sec  of  equal  ly  -  spaced  points  in  coordinate  space,  and  the 

(2) 

Laplacian  terra  for  fig^  can  also  be  expressed  by  finite  differences.  Then  we  get 

a  set  of  inhomogeneous  linear  equations,  which  can  be  written  in  matrix  form  and 

diagonalized  by  the  standard  Gauss  elimination  technique.  By  making  a  good  choice 

(2) 

for  the  initial  g^  (r)  .  the  rate  of  convergence  will  be  rapid,  so  that  the  optimized 
distribution  function  can  be  readily  obtained. 


III.  Numerical  Results  and  Discussion 

In  this  section  we  apply  the  results  obtained  in  previous  sections  to  tvo- 

3  3  4  3 

dimensional  He  and  (Bose)  He-  He  mixtures  in  the  zero  He  concentration  limit.  To 


be  realistic  in  the  two-dimensional  mixture,  It  is  necessary  to  consider  the  He 

4  19  3 

coverage  of  a  few  atomic  layers  on  the  He  background.  At  T  ■  0  K  a  He  particle 

,  ,  4 

behaves  nearly  like  He  particle,  and  thus  wc  assume  that  the  potential  V  (r) 

QP 


11 


appearing  in  the  Hamiltonian  [ Eq .  (2.1)]  are  all  of  the  form  of  a  Lennard- Jones 

potential  given  as 


Va.(r) 


M(;)12  *  (J)6] 


(3.1) 


(2 ) 

with  e  -  10.22  K  and  a  -  2.556  A.  As  an  initial  trial  guess  for  g22  (r) .  we  adopt 
(2) 

the  results  for  g^2  (r)  and  S22(k)  obtained  from  the  Jastrow  wavefunction,  which  are 
only  known  for  r  <  5  A  and  k  <  5  A  \  such  that  the  short-range  behavior  of  g(r)  at  a 
density  of  0.03  A  coincides  with  the  results  obtained  from  the  PPA.  The  initial 
form  of  g22^ (t)  is  evaluated  using  an  iterative  scheme,^  since  the  long-range 
behavior  of  g(r)  depends  on  the  value  of  S(k)  for  small  k,  and  the  exact  value  of 

S(k)  for  k  >  5A  ^  depends  on  the  value  of  g(r)  for  small  r.  The  initial  value  of 

(2)  (2) 

g^2  (r)  is  determined  by  optimizing  g^2  (r). 

(3) 

For  the  triplet  correlation  functions  we  assume  that  f  (r)  has  the  same  form 
for  all  the  triplets  and  choose  the  parametrized  form  given  as 


r(3) 


f  (ri,r2’r3>  ”  GXpt  /  A^(rij)C<rik^r  r 


ij 


'rik] 


cyclic 


(3.2) 


with 


f (r)  -  F(r)exp{[(r-rt)/w]2)  (3.3) 

F(r  <  b)  -  [ (r  -  rb)/rb]3  ,  F(r  >  rfa)  -  0  . 

*  . 

Here  r..  and  r..  are  unit  vectors  in  the  direction  of  radial  vectors  r.,  and  r.,  , 
lj  lk  lj  ik’ 

respectively.  We  choose  the  parameters  rb>  rt,  A  and  u  as  those  which  are  adopted  in 


12 


the  three-dimensional  case,  and  note  that  the  energies  and  radial  distribution 
functions  are  not  significantly  affected  by  the  choice  of  the  above  parameters. 
Taking  the  terms  of  elementary  diagrams  up  to  fifth  order  with  the  triplet 
correlation  factor  and  approximating  E^Cr)  by  E^Cr),  the  calculation  of  the 
in  t e gr o  -  di  f  f  e r ent ial  equation  {Eq.  (2.24)]  is  carried  out  by  using  the  Gauss 
elimination  technique  until  r  -  20  A. 

The  ground-state  energies  per  particle  for  pure  liquid  ^He  and  ^He-^He  mixtures 

are  listed  in  Table  1.  By  considering  the  triplet  correlation  factor  and  elementary 

diagrams,  the  energies  are  lowered  by  about  21-96%  compared  with  those  considering 

the  two-body  factor  only,  and  the  equilibrium  density  increases  for  the  pure  case  but 

decreases  for  the  mixture  case.  The  changes  by  the  triplet  correlation  and 

elementary  diagrams  are  shown  in  Figure  1.  For  the  pure  case  (solid  line),  the 

.  -2 

ground-state  energy  is  -0.646  K  at  the  equilibrium  density  of  p  -  0.0350  A  and 

18  20 

similar  to  the  results  of  the  BBGKY,  the  molecular  dynamics  and  the  Monte  Carlo 

4  .  „  o 

studies.  For  the  mixture  case  (dashed  line),  the  equilibrium  density  is  0.0336  A 

with  a  ground-state  energy  -0.542  K.  The  rate  of  drifts  in  energies  are  similar  to 

that  in  the  three-dimensional  system.  Moreover,  we  can  view  the  chemical  potential 
3  U 

of  He  at  T  =  0  K  as  independent  of  the  density  of  the  He  background  due  to  complete 
phase  separation.^ 

We  note  that  one  of  us  (CIU)  has  evaluated  the  ground-state  energies  of  one-  and 

21  99 

two-dimensional  Bose  liquids  and  charged  boson  systems  using  ring  diagrams  c  and 

23 

the  self-consistent  field  approximation.  In  these  results  the  ground-state 
energies  of  one-  and  two-dimensional  Bose  fluids  consist  of  repulsive  and  attractive 
parts  with  the  adoption  of  a  soft  potential  with  a  Lennard- Jones  tail.  For  our 
current  two-dimensional  boson  system,  the  ground-state  energy  is  proportional  to 
-  2/3 

R  with  an  additive  positive  constant,  where  R  is  the  ratio  of  the  mean 

cl  cl 
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separation  between  particles  to  the  Bohr  radius  within  the  self-consistent 

approximation.  This  result  is  different  from  that  of  the  ring  diagram  approximation, 

which  does  not  include  an  additive  positive  constant. 

Figure  2  illustrates  the  optimal  radial  distributions  of  g^9  (r)  (dashed  line) 

(2) 

and  g12  (r)  (  solid  line)  at  the  respective  equilibrium  density,  while  their  numerical 

values  at  various  density  are  listed  in  Table  2.  As  the  density  increases,  the 

(2)  (2) 

maximum  values  of  (r)  and  §2.2  increase.  At  the  equilibrium  density  the 

(2) 

maximum  position  of  g^  (r)  shifts  to  the  right  (direction  of  increasing  radial 

(2) 

distance)  in  comparison  with  the  fixed  g^  (r)-  This  behavior  is  due  to  the  fact 

3  A 

that  the  volume  swept  by  He  is  larger  than  that  by  He.  The  maximum  value  of 
( 2 )  (2) 

g^2  (t)  is  smaller  than  that  of  (r)  and  broadened  around  the  maximum  point 

3 

because  of  less  localization  of  the  He  impurity  with  a  larger  zero-point  energy.  In 

( 2 )  23 

this  sense  the  behavior  of  g^^' (r)  agrees  with  the  three-dimensional  case.  With 

the  consideration  of  the  triplet  correlation  function,  the  maxima  of  g^iO  and 

(2) 

g22  (r)  increase  but  are  hardly  affected  by  elementary  diagrams.  However,  the 

ground-state  energies  are  lowered  not  only  by  triplet  factors  but  by  elementary 

21  22  25 

diagrams.  We  note  that  the  optimized  radial  distribution  functions  ’  '  of  one-, 

2  ^  A 

two-  and  three-dimensional  hard-sphere  Bose  systems  decrease  as  r  ,  r  and  r* 
respectively,  at  short  distances,  while  at  long  distances  they  oscillate  and  tend  to 
go  to  unity  as  r  -»  ®,  which  is  due  to  the  contribution  from  the  core  part  of  the  soft 
potential.  These  three  proportionalities  represent  the  main  long-distance  behaviors 
of  phonons  in  each  dimensional  Bose  liquid. 

Figure  3  shows  the  optimized  liquid  structure  functions  at  various  densities. 
The  initial  slope  decreases,  which  implies  a  decrease  of  sound  velocity,  but  the 
maximum  increases  and  shifts  to  the  right  with  increasing  density.  A  comparison 
between  S^Ck)  (solid  line)  and  S22^k^  *  1  (dashed  line)  at  equilibrium  density  is 

3 

shown  in  Figure  4.  The  addition  of  the  impurity  He  leads  to  a  decrease  of  the 
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maximum  and  shift  of  the  maximum  position  to  the  right  in  k- space.  This  behavior 

obtained  is  as  expected  in  a  scattering  problem  and  agrees  with  that  of  the  three- 
25 

dimensional  case  except  in  the  region  of  k  -*  0. 

In  Figure  5  we  plot  the  bare  Lennard- Jones  potential  V  the  induced  potential 

W12^t0  (V.^)  anci  effective  potential  According  to  Table  3,  we  observe  for 

the  effective  potential  that  as  the  density  increases,  its  hard  core  radius  (V  - 

ef  f 

0)  decreases  and  potential  becomes  deeper,  while  the  repulsive  barrier  height  grows 
at  about  the  same  position. 

IV.  Conclusions 

3  A 

We  have  investigated  the  ground-state  properties  and  behaviors  of  the  He-  He 

system.  The  introduction  of  the  triplet  correlation  and  elementary  diagrams  up  to 

fifth  order  in  the  HNC  approximation  leads  to  successful  improvement  in  the  ground- 

state  energy  and  optimized  two-body  distribution  functions.  Furthermore,  the  triplet 

correlation  factor  contributes  to  the  change  of  both  energy  and  the  distribution 

function,  while  elementary  diagrams  do  not  change  the  distribution  function.  By 

3  4 

addition  of  an  impurity  He  into  a  He  background,  the  energy  increases,  the 

distribution  function  moves  to  the  right,  and  its  maximum  decreases  and  is  broadened 

around  the  maximum  position,  due  to  larger  zero-point  energy  and  less  localization  of 
3 

the  He  particle. 
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Table  1.  Ground-state  energies  due  to  the  two-body  factor  only  (2) ,  two-  and 
three-body  factors  (2  +  3)  and  two-  and  three-body  factors  plus 
elementary  diagrams  (tot)  for  the  pure  and  mixture  cases  at  various 
densities . 


p 

(A'2) 

E0.n,l*<2> 

00 

E_  .  (2+3) 

0  ,mix 

(K) 

En  .  (tot) 

0 ,  mix 

(K) 

E0.pure<2) 

(K) 

E„  (2+3) 

0 ,  pure 

00 

E_  (tot) 

0 .pure 

(K) 

0.0250 

-0.28 

-0.38 

-0.45 

-0.39 

-0.45 

-0.50 

0.0275 

-0.30 

-0.40 

-0.48 

-0.44 

-0.48 

-0.55 

0.0300 

-0.32 

-0.43 

-0.50 

-0.48 

-0.53 

-0.58 

0.0325 

-0.35 

-0.45 

-0.53 

-0.46 

-0.56 

-0.62 

0.0350 

-0.38 

-0.47 

-0.54 

-0.46 

-0.59 

-0.65 

0.0375 

-0.35 

-0.45 

-0,52 

-0.43 

-0.58 

-0.63 

0.0400 

-0.33 

-0.44 

-0.51 

-0.40 

-0.57 

-0.61 
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Table  2.  Radial  distribution  functions  at  various  densities. 


P  -  o. 

, 02 j  A'2 

P  -  o. 

.030  A'2 

p  -  0. 

.035  A'2 

r  (A) 

g22  ^  ^ ^ 

g12(r) 

g22(r) 

Sl2(r> 

g22 

s12(r) 

2.0 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

2.2 

0.010 

0.010 

0.010 

0.011 

0.010 

C.011 

2.4 

0.067 

0.065 

0.070 

0.066 

0.071 

0.069 

2.6 

0.208 

0.207 

0.211 

0.209 

0.215 

0.209 

2.8 

0.409 

0.404 

0.418 

0.408 

0.425 

0.413 

3.0 

0.628 

0.608 

0.640 

0.621 

0.656 

0.630 

3.2 

0.835 

0.801 

0.839 

0.815 

0.841 

0.828 

3.4 

0.988 

0.960 

0.992 

0.963 

0.997 

0.970 

3.6 

1.086 

1.055 

1.095 

1.066 

1.101 

1.075 

3.8 

1.140 

1.109 

1.159 

1.117 

1.160 

1.126 

4.0 

1.188 

1.140 

1.189 

1.155 

1.198 

1.160 

4.2 

1.181 

1.145 

1.182 

1.158 

1.190 

1.162 

4.4 

1.169 

1.131 

1.168 

1.135 

1.173 

1.140 

4.6 

1.121 

1.095 

1.121 

1.099 

1.123 

1.103 

«*.8 

1.084 

1.077 

1.084 

1.078 

1.085 

1.075 

5.0 

1.050 

1.048 

1.043 

1.043 

1.043 

1.040 

5.2 

1.015 

1.019 

1.015 

1.018 

1.010 

1.011 

5.4 

0.990 

0.996 

0.987 

0.992 

0.986 

0.990 

5.6 

0.972 

0.978 

0.969 

0.974 

0.967 

0.972 

5.8 

0.957 

0.963 

0.956 

0.961 

0.955 

0.959 

6.0 

0.951 

0.956 

0.948 

0.953 

0.947 

0.952 
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Table  3.  Induced  potentials  ^  and  effective  potentials  at 


various  densities, 


p  -  0. 

025  A*2 

P  -  0. 

030  A'2 

p  -  0.35 

A’2 

r(A) 

Vind<K> 

Veff(K> 

V.  J 
ind 

Veff 

V.  ^ 
md 

Veff 

2.4 

-1.056 

26.332 

-1.625 

25.763 

-2.401 

24.987 

2.6 

-0.756 

-4.347 

-1.162 

-4.754 

-1.752 

-5.344 

2.8 

-0.480 

-10.447 

-0.942 

-10.709 

-1.159 

-11.126 

3.0 

-0.226 

-9.882 

-0.367 

-10.023 

-0.623 

-10.279 

3.2 

0.009 

-7.850 

-0.027 

-7.887 

-0.141 

-8.000 

3.4 

0.231 

-5.816 

0.290 

-5.758 

0.302 

-5.744 

3.6 

0.437 

-4.129 

0.579 

-3.987 

0.702 

-3.863 

3.8 

0.607 

-2.828 

0.812 

-2.624 

1.024 

-2.412 

4.0 

0.757 

-1.828 

1.026 

-1.568 

1.315 

-1.279 

4.2 

0.950 

-1.021 

1.279 

-0.692 

1.643 

-0.329 

4.4 

1.216 

-0.294 

1.637 

0.126 

2.086 

0.575 

4 . 6 

1.546 

0.378 

2.052 

0.884 

2.588 

1.421 

4.8 

1.828 

0.917 

2.364 

1.453 

2.964 

2.053 

5.0 

1.885 

1.168 

2.359 

1.642 

2.960 

2.244 

5.2 

1.625 

1.057 

1.968 

1.490 

2.501 

1.932 

5.4 

1.123 

0.669 

1.328 

0.874 

1.745 

1.291 

5.6 

0.576 

0.210 

0.694 

0.328 

0.995 

0.629 

5.8 

0.154 

-0.143 

0.249 

-0.048 

0.466 

0.168 

6.0 

-0.088 

-0.331 

-0.022 

-0.221 

0.191 

-0.052 
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Figure  Captions 

1.  Ground-state  energies  of  the  pure  ^He  system  (solid  line)  and  (Bose)  He-  He 
mixture  in  the  impurity  He  limit  (dashed  line)  as  a  function  of  density.  E(2) , 
E(2+3)  and  E(tot)  represent  the  energy  containing  the  two-body  correlation  factor 
only,  energy  by  two-  and  three-body  factors,  and  total  energy  by  the  above 
factors  plus  elementary  diagrams,  respectively. 

2.  Optimized  radial  distribution  functions  g£^(r)  (solid  line)  and  (r)  (dashed 

line)  at  equilibrium  density. 

3.  Optimized  liquid  structure  function  S12(k)  at  various  densities. 

4.  Optimized  liquid  structure  function  S^2(k)  (solid  line)  and  S22(IC)"1  (dashed 
line)  at  the  equilibrium  density. 

Bare  Lennard- Jones  potential  (V^j),  induced  potential  (v^ncj)  anc*  effective 
potential  (ve£f)  aC  the  equilibrium  density. 


5. 


